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ABSTRACT 


The stability of an incompressible, inviscid, electrically conducting fluid cylinder in the presence 
of axial current is considered taking into account the effect of surface tension. The effect of gravity 
is taken into account approximately. The results are found to be in reasonable agreement with 
experiment. 


1. Introduction 


In recent years the theory of stability of constricted gas discharges was studied in 
great detail in view of its importance in the problem of controlled production of 
thermonuclear energy. Kruskal and Schwarzchild [1], Taylor [2, 3], among others, 
showed that a pinched gas discharge is highly unstable for “‘sausage’’ and ‘“‘kink”’ 
deformations. Experimental evidence [4] confirms that theoretical results are quali- 
tatively correct. 

In view of the importance of the stability problem, one would like to have quanti- 
tative comparison between theory and experiment under suitable conditions. With 
this end in view a liquid scale model of pinched gas discharge has been constructed by 
Dattner, Lehnert and Lundquist [5]. In their experiment a mercury jet of 14 cm 
length and 0.2 cm radius, falling vertically from a circular hole in a tank, plays the 
role of a pinched gas discharge. It is observed that when an axial current is passed 
through the mercury jet without any superimposed axial magnetic field, the jet 
develops a sausage instability while in the presence of axial magnetic field and axial 
current it develops corkscrew instability. This experiment is not only free from the 
difficulties of pinch discharges, but also has the merit that it allows quantitative 
comparison with theory while bringing forth the qualitative information about insta- 
bility of pinched gas discharges. It is the aim of this paper to develop a theory for 
axi-symmetric perturbation in the presence of an axial current without an axial 
magnetic field and compare the results of the theory with the above-mentioned 
experiment. The effect of an axial external magnetic field will be considered in a 


later paper. 
2. Basic assumptions and equations 


Let us consider (as a working model) an infinitely long fluid cylinder of circular 
cross-section of radius R, surrounded by vacuum. We assume that the fluid is in- 


1 On deputation from The Tata Institute of Fundamental Research, Apollo Pier Road, 
Bombay, India. 


16:3 241 


c. s. MURTY, Instability of conducting fluid cylinder 


compressible and inviscid and has a finite electrical conductivity. A current J is 
assumed to be flowing along the axis of the cylinder. In the following analysis we 
shall first ignore the effect of gravity, but discuss this effect at the end. 

In C.GS. electromagnetic units (ignoring the effect of displacement current) the 
basic equations are 


9 = -vp+ ix, (2.1) 

V:v=0, } (2.2) 

VxH=4nj, (2.3) 

ae ~yxE, (2.4) 

j—o(E+v~xH), (2.5) 

ee Hee 0), (2.6) 

and I =constant. (2.7) 


The symbols in the eqns. (2.1) to (2.6) have their usual meaning [6]. In the last eqn. 
(2.7) J =total current passing along the fluid cylinder. All the eqns. (2.1) to (2.6) are 
the basic equations of any problem in hydromagnetics of an incompressible fluid of 
finite electrical conductivity. Eqn. (2.7) is a constraint on the system in the present 
case denoting that the source of current is outside fluid cylinder and supplies a current 
of constant value. 


Equilibrium 
The equilibrium state of the fluid cylinder can be described by 
0=-Vpet+jexH., (2.8) 


and eqns. (2.3) and (2.7). If jp is the equilibrium current density along the axis of the 
fluid cylinder which we assume coincides with z-axis of circular cylindrical coordi- 
nates, we have 


Jo= Re: (2.9) 


One can show that the pressure, p;, in the interior of the fluid cylinder at a distance 
r < R, from the axis is given by 


—z < 2 / p22 
ss ahie  Picimls 4 (Ro—r*), (2.10) 


where T' = surface tension of the fluid (mercury in our problem). 
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3. Criterion of instability of the equilibrium for axi-symmetric perturbation and 
the method of solution 


In order to determine whether the equilibrium state given by eqns. (2.9) and (2.10) 
is stable or not, we assume that the fluid cylinder is subjected to an axi-symmetric 
perturbation which deforms the surface of the cylinder into a shape given (in circular 
cylindrical co-ordinates with z-axis coinciding with the axis of the cylinder) by 


r= K+ a cos kz, (3:1) 
where a is a function of time subject to the condition 
a<k, (3.2) 


so that we can neglect terms proportional to a?. Due to the fact the fluid is assumed 
to be incompressible we find that 
2 


R= +S, se ae hy (33 


so that one can ignore the difference between R, and R. If we assume further that @ 
is periodic in time so that 
a =a, cos wt, (3.4) 


the problem of investigating the stability of the fluid cylinder will be shifted to finding 
the values of w. If w is found to have an imaginary value for any real value of k 
we conclude that the fluid cylinder is unstable for perturbations that deform its sur- 
face into a shape given by eqn. (3.1). The values of w can be found from eqn. (2.1) 
if we know the distribution of pressure, p, as well as current density and magnetic 
field in the perturbed condition. These quantities will be calculated in the following 
sections. 


4, Distribution of perturbed current density and magnetic field 


The fundamental quantity in the calculation of w is the perturbed current density 
distribution. We shall calculate this quantity in the “magnetostatic approximation’. 
This approximation is good when the electrical conductivity of the fluid is so small 
that the characteristic time for diffusion of magnetic lines of force is smaller than the 
characteristic time for fluid motions [7]. A more direct way of looking at this approxi- 
mation is as follows. 

If we combine eqns. (2.4) and (2.5) we obtain 


of = — vx [Lv], (4.1) 
o 

If the electrical conductivity o of the fluid is infinitely large eqn. (4.1) leads to the 
well-known equation of induction for frozen magnetic lines of force: 


oH 


au Vv x(vxH). (4.2) 


243 


c. s. MuRTY, Instability of conducting fluid cylinder 


If, on the contrary, the electrical conductivity of the fluid is vanishingly small we 


find that 
Vy xj=0, (4.3) 


so that the current density can be derived as a gradient of scalar function, say, p: 
j=-Vy- (4.4) 


As far as magnetohydrodynamic problems are concerned mercury can be considered 
as a poor electrical conductor and therefore we assume that eqn. (4.4) is valid. By 
virtue of eqn. (2.3), y satisfies the Laplace equation: 


vV2y =0. (4.5) 


In order to determine uniquely the perturbed current density distribution, we 
impose the constraint given by eqn. (2.7) in an equivalent form. We demand that 
there is no current flowing along the normal to the perturbed boundary, or 


(j X D)r=R+ac0s ke = 9, (4.6) 


where n is the unit normal to the surface of the deformed cylinder. One can show that 


n= (1, 0, 0) + (0, 0, ak) sin kz. (4.7) 


We remark here that eqns. (4.3) to (4.6) are formally the same as those considered 
by Chandrasekhar and Fermi [8] in connection with stability of an infinitely conduct- 
ing incompressible fluid cylinder with an axial magnetic field frozen in the medium. 
Following their method one can show that the perturbed current density in the 
present problem is given by 


1 =(j+,.9, Je), ] 
: aed f 
where jr = —akj, 7 v9 sin kz, 
1 (k &) (4.8) 
STS inlo(ber 
and Is = jo 7 2KIo ae cos kz, 
» | 


where I,,(kr) is the modified Bessel function of order n. Using eqns. (2.3) and (4.8) we 
find that the perturbed magnetic field distribution is given by 


H= (0, Hp, 0), (4.9) 
where 
y ~Ly(kr) 
2 =A u 
Hpelig en 


a= 


coskz, for r<R+acos kz, 
(4.10) 


2ajoR?/r, for r>R+acoskz. 
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5. Distribution of pressure inside the deformed cylinder 


The knowledge of the distribution of current density inside the deformed cylinder 
enables us to calculate the distribution of pressure. If we take the divergence of both 
sides of eqn. (2.1) and use eqns. (2.2), (2.3), (4.4) and (4.8) we find, retaining terms 
up to first order, that 


I, (kr) 


2 = —4 2 2 ‘ Zz. 5. 
VP JO ce RU RTT cos k (5.1) 
In order to solve eqn. (5.1) let us write 
p=, (r)+ap, (r) cos kz. (5.2) 


If we substitute this expression for p in eqn. (5.1) and separate terms of corresponding 
orders we obtain two inhomogeneous differential equations which can be solved by 
standard methods. We only give the result which shows that 


Po (r)=A-axjor, 
I,(kr) p> (5.3) 
I, (kR) 


and P(r) =BIy (kr) + 4a jer 


where A and B are arbitrary constants to be determined by boundary conditions. 

To determine the constants A and B we note that the pressure, p,, at a point just 
inside the deformed boundary depends only on the surface tension and curvature 
of the surface since we assumed that the pressure outside the fluid cylinder is zero. 
One can show that [9] 


(k? R? —1) cos kz. (5.4) 


1p EY 7 


The boundary condition we impose on the solution of eqn. (5.1) is that the right-hand 
side of eqn. (5.2) evaluated at r = R +a cos kz equals the right-hand side of eqn. (5.4). 
We then find that 


Tes <3 
= R 
A RI ’ 
(5.5) 


Tete et, 
= a3 (R= 1)=2 aR | 
a 2 Tem [we RY 2a8 


Therefore, eqns. (5.2), (5.3) and (5.5) define uniquely (apart from a trivial additive 
constant) the distribution of pressure inside the deformed cylinder in the range 
0<r<R-+<a cos kz up to terms of order a. 


6. Dispersion relation 


The calculation of @ is relatively simple matter since we obtained all the quantities 
on the right-hand side of eqn. (2.1). If we use eqn. (3.4) we can write (neglecting terms 
of higher orders) the radial part of eqn. (2.1) as 
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Op: | 
—eo° ={—2-+(jx HB), , 6.1 
ow acos kz ap ae (6.1) 


where the suffix to the brace on the right-hand side of eqn. (6.1) indicates that the 
quantities inside that brace must be evaluated at the deformed boundary. Using 
eqns. (4.8), (4.10), (5.2), (5.3) and (5.5) we find that eqn. (6.1) gives the required 
dispersion relation: 


are nfiel-axfies eB, ag 

where x=kR, 
N= /nRT, (6.3) 
and OQ? = ow" R°/T. (6.4) 


We can draw the following conclusions from eqn. (6.2). 

(i) Q2 is a real number for real values of & which means that there are no states of 
overstability, a result proved by Bernstein et al. [10] in the case of infinitely conduct- 
ing plasma. 

(it) By taking the asymptotic values of the Bessel functions one can show that 


Lt Q?—> — 2? (1+ N)/2, 


z30 


and Lt Q?> + 0, 


t—>0o 


for all finite values of N. Since the modified Bessel functions are monotonically 
changing functions of the argument we can conclude from the above asymptotic 
limits that Q? has one more zero, apart from the one it has at the origin, at some value 
of x which we may denote by ~,. In the range 0 < x < a,, we find that Q? < 0 and has a 
minimum at some value of x which we denote by z,,. Therefore in the range 0 < x < 2, 
the values of m are pure imaginary and hence we conclude that for disturbances that 
impress on the cylinder wavelengths corresponding to the range 0 < x < 2,, the fluid 
cylinder is unstable. For values of wavelength corresponding to # = z,, the instability 


Table 1. Variation of x,=27R/A, and %,=227R/A_ with N = I?/x RT. 


: 
The values of x, and a,, are accurate up to +0.02 due to graphical method of solution adapted 
for obtaining them. 


N | iy oe | Le 
0 1.00 0.695 1.44 
1 1.37 0.94 1.46 
2 1.62 1.10 1.47 
3 1.83 1.24 1.48 
4 2.01 1.34 1.50 
5 el Ne 1.44 L51 


a ee ee 8 Med A el nt 
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° Experimental values 


Fig. 1. Variation of growth constant 
| @,,|(ms)-! with axial current I (Amp) 
for mercury jet of 0.2 cm radius; 7 = 487 
dyne cm—'; @ = 13.6 g cm-*. Axi-sym- 
metric perturbation. 


grows at a maximum rate. We further note that the term inside the brace of eqn. 
(6.2) is always positive, thereby contributing a negative quantity to Q? since NV > 0. 
This shows the unstabilising effect of current flowing along the axis of the fluid 
cylinder. 

For a given value of N the values of x, and w,, are obtained graphically and the 
results are shown in the Table 1. 

Corresponding to the value of x = x,, the value of |«,,| is calculated as a function 
of current J (amp) for mercury jet of 0.2 cm radius (7' = 487 dyne em-! and 9 = 13.6 g 
em-*) and the results are shown in Fig. 1. In the same figure the values of |@| 
measured from the film taken by Dattner [5] are also plotted. The experimental 
values are consistently larger than the theoretical values. 

There are mainly three possible reasons for the discrepancy between the experi- 
mental and theoretical results. 

(1) The actual surface tension of mercury can be smaller than the value used in 
the theoretical calculations on account of surface impurities. 

(2) The effect of viscous forces is neglected in theory. This is a good approximation 
in the present problem since the effect of viscous forces is several orders of magnitude 
smaller than electromagnetic forces in the case of mercury [11]. 

(3) The jet is falling vertically under gravity. This is expected to appreciably affect 
the stability and therefore will be discussed in some detail in the next section. 


7. Effect of gravity 


It has been assumed in the theory developed so far that the fluid is unaffected by 
gravity. If the fluid is falling vertically downwards due to gravity we have to take 
into account (1) the decrease of cross-section of the fluid cylinder and (2) the lengthen- 
ing of undulations that develop due to instability as the velocity increases. As the 
jet of fluid emerges from the orifice the disturbance that impresses upon it undula- 
tions of length approximately equal to 0.8 times the circumference of orifice (when J 
= 300 A) are most unstable. But as the jet moves faster the parameter «= 22 R// 
decreases due to effects (1) and (2) and therefore the value of m changes from its 
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maximum value near the orifice.! The effect on w due to variation of radius is more 
prominent at initial stages near the orifice and for smaller initial velocities of the jet. 
At large distances from the orifice the effect on w due to the variation of radius is 
relatively smaller than that due to variation of wavelength. 

Due to the lengthening of undulations, the curvature of the surface of the jet also 
changes. We shall ignore this effect, but in a complete theory one should include this 
also. 

It is possible to describe approximately the effect of gravity by considering the 
variation of wavelength 4 and mean radius R along the axis of the jet, and defining 
an “instantaneous growth constant”’ @inst at a given distance / from the orifice. If we 
ignore the variation of R in distances less than or equal to one wavelength, one can 
show (for distances | >A from orifice) that wavelength /(/) is given by 


AD=—49T5+ [29 TO (h+))}, (7.1) 


and the radius R(l) is given by 


1]4 
RQ=R,/|1 +5 ; (7.2) 


In eqn. (7.1) 7’, is defined as the time interval in which the fluid emerging from the 
orifice moves a distance equal to /,, or 


Am =T,(2gh)' +4973, (7.3) 


where h = the height of the fluid above orifice and g = the acceleration due to gravity. 
T, is calculated from the theoretical value of 4,, defined in the gravity free case. In 
the experiment of Dattner [5] the value of h is about 3 cm. Taking h = 3 cm the value 
of Wins, for various values of | is calculated for J =300 A and the results are shown 
in Table 2. The large initial value of wns_ is due to the rapid change of radius of the 
jet as it emerges from the orifice. The decrease of @nst at larger distances from orifice 
is due to the increase of wavelength which dominates the effects due to variation in 
radius as may be seen from eqns. (7.1) and (7.2). 

Calculations show that the above results are qualitatively the same when the cur- 
rent is zero. This should be contrasted with the conclusion of Rayleigh [12], namely, 
that “‘the change of conditions during fall (under gravity) has a protective influence’. 
This is true only at great distances from the orifice but at closer distances the changing 
conditions have an unstabilising effect. 


Table 2. Variation of ‘instantaneous growth constant”, winst, due to gravity along 
the mercury jet. 


1 is the distance from the orifice of radius 0.2 em, the axial current is 300 A. 


1 (cm) | 2 | 4 | 6 | 8 


inst * 10° (ms)* 


1 In this section we shall write @ in place of |@|. 
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It is interesting to note that the “instantaneous” value of w at 1 between 6 cm and 
12 em varies only by about 8% with a mean ~ 0.088 (ms)-1 (refer Table 2) which 
should be compared with experimental value of 0.099 + 0.010 (ms)— for J = 300 A. 
(Fig. 1). We remark here that the experimental value is the mean result of observa- 
tions made in the lower half of the jet (total length 14 cm). 

Finally, we point out that the value of «,, at the orifice of the jet can be determined 
experimentally by observing the variation of wavelength along the fluid cylinder and 
extrapolating to the orifice. As the deformation becomes observable only in the lower 
half of the jet the extrapolation is a bit uncertain but nevertheless it is found to give 
reasonable agreement with theoretical value. 


8. Concluding remarks 


We have shown that a current passing through a conducting fluid cylinder has an 
unstabilising effect. In the absence of surface tension the fluid cylinder will break 
into pieces of very small length which should be compared with the instability of 
pinch discharge. The force due to surface tension has a stabilising effect at small 
wavelengths. When the fluid is falling under gravity, neglecting the effect of the 
variation of curvature of the deformed surface, the effect of variation of wavelength 
and radius can be approximately described by eqns. (7.1) and (7.2). Viscosity has 
been neglected but is expected not to change appreciably the above results in the 
case of mercury jet. 
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